DENSITY FLUCTUATIONS IN THE GALACTIC HALO 
AND EXPERIMENTAL SEARCHES FOR DARK MATTER 
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Clumping of elementary dark matter in the Galaxy halo may be inevitable. If so, the nondetection of certain 
dark matter candidates could simply mean that the local halo density is low. Conversely, indirect detection 
] of annihilation products could be facilitated, perhaps to an embarrassing degree. 
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' The interpretation of the results of experiments designed to search for elementary particle dark matter in 

the Milky Way halo relies on specific suppositions concerning the local mass density and velocity distribution 
j^jrj' of the unseen particles. Predicted detection rates 1 have all been computed under the assumptions that the 

halo density profile is smooth, and the velocity distribution is at least approximately Maxwellian, with a 
one-dimensional velocity dispersion a = V c /V% ~ 160 kms -1 , where V c ~ 220 kms -1 is the Galactic rotation 
speed. However, an experiment that operates for a time T ~ 1 yr only probes a region that extends to a 
distance D ~ aT ~ 5 x 10 14 cm 35 A. U., which is about the distance to the planet Pluto; for comparison, 
the distance to the center of the Galaxy is Dq « 8.5 kpc « 2 x 10 9 A.U. While few would doubt that the 
■ halo density may be regarded as smooth on some intermediate length scale, say Ss 0.01 — 1 kpc, there is 

reason to question whether the halo really ought to be perfectly smooth on length scales of order the size of 
the inner Solar System, ~ 10~ 8 — 10~ 6 kpc. 

If the halo is clumpy on small scales, then it is possible that presently the Earth is not immersed in a 
bath of dark matter particles, so that the non-detection of such particles need not imply their absence. As 
\Q , a simple but extreme example, suppose that the dark matter particles are organized into lumps of internal 

ON ■ density p c i 3> ~p, where p w O.OlAf© pc -3 is the mean halo density locally; then the probability that the 

Earth sits in a dark matter clump is the volume filling factor, / = ~p/p c i *C 1. If instead the halo contains 
O ^ a continuous mass spectrum of dark matter lumps, so that ~pip(M)dM is the contribution to the local halo 

density from clumps of mass M and internal density p c i(M), then the volume filling factor of the halo is 

H 

f = JS I dM4>(M)pJ(M), (1) 

which may also be very small if the dark matter is concentrated into dense clumps. Consequently, experimen- 
talists searching for elementary particle dark matter in a clumpy halo could simply be unlucky. Conversely, 
a fortunate experimenter bathed in a dark matter cloud could detect halo particles at an incredibly large 
rate compared to predictions based on the assumption of a smooth halo. 

Although the problem of small scale halo structure is generic, we focus here on a specific case: neutralinos 
of mass m < GeV whose annihilation is predominantly s- wave. 2 The cosmology of such particles is governed 
by m, and a dimensionless parameter, C = Mpim{<jv), where Mpi is the Planck mass and (av) is the 
annihilation rate coefficient for the fermions. Annihilation equilibrium fails at a temperature Tf ~ m/lnC, 
allowing a significant present day abundance, amounting to a total mass density po ~ mTJ 3 In C /C, where 
To is an appropriate present day temperature (e.g. about 2°K for neutrinos); at density parameter f2o and 
Hubble constant H = 100/iokms" 1 Mpc -1 , C ~ 10 9 m(GeV)(n hl)' 1 . The neutralinos remain in thermal 
contact with the ambient gas of annihilation products until the temperature falls to Tk ~ toC" 1 / 4 , after 
which they decouple completely from the rest of the cosmological soup. 3 

After Tk, the neutralinos stream freely along straight line paths. Consequently, any fluctuations 
in the density of dark matter particles will be erased on scales smaller than the free-streaming length, 
which corresponds to a mass 3 M fs ~ {M Pl /m 2 ) In C/C 5 / 8 - lO- 4 M (fi o ^o) 5/8 /[ m ( GeV )] 21/8 - More- 
over, by Liouville's theorem, the microscopic phase space density is preserved after Tk, when its value 
was A p ~ In C/C 11 / 8 ~ 10 _11 [rio/io/m(GeV)] 11 / 8 . The coarse-grained phase space density of the halo is 
A/j ~ 10 _32 /[m(GeV)] 4 <§; A p , which suggests that some fine-grained clustering of neutralinos in the halo is 
inevitable. 
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The nature of the clustering is still not certain. One possibility is that dark matter clumps of mass Mf s , 
which are expected to be among the first objects to condense cosmologically, survive intact to the present 
day and dominate the Galactic halo. Research along these general lines has been reported at this conference 
by Gurevich and collaborators, who also suggest that non-baryonic clumps may be the MACHOs. Another 
possibility is that clumps are not permanent but come and go in a dynamical fashion as the result of small 
scale 'gravitational turbulence' in the Galaxy halo: locally space-time dependent density fluctuations with 
possibly stationary statistical properties. 

A simple picture that could lead to a stationary statistical state for the halo would involve binary 
fragmentation and merger of dark matter clumps in the halo. To be concrete, let us focus on a discrete 
model in which Nk is the number of clumps of mass Mk = Mo2~ k , where k is an integer; in a continuous 
model, Nk would correspond to tp(M) in Eq. (1). Suppose that fragmentation of a cloud of mass Mk results 
in the production of a pair of bound clouds of mass Mk+i; merger of a pair of clouds with equal masses Mk 
forms a cloud of mass Mk-i- The net rate of formation of clouds of mass Mk is then 

^ = -{v k + 2v k -i)Nk + 2v k -iNk-i + VkNk+i, (2) 

where Vk is the rate of fragmentation of a mass Mk and Vk is the rate of merger to form a mass Mk- If we 
assume that Vk ~ \[Gp~k and Vk — cvk, where pk is the internal density for mass Mk and c is independent 
of scale, then Eq. (2) has the steady state solution Nk oc M^° S2C ; notice that this is independent of the 
dependence of density on mass scale, and that the mass spectrum is tilted toward small masses for c < 1/2. 
(This steady-state solution only requires the assumption that v c = cvk, not the additional physical conjecture 
that the characteristic merger and collapse rates are approximately free-fall.) Assuming that the basic physics 
of collapse and merger is scale free also leads to pk oc M^ 2ri ; if fragmentation tends to leave behind bound 
pairs, then 77 < 1 in general. 4 The mass spectrum extends down to the 'Jeans mass', where the internal phase 
space density of a lump is ~ A p ; this is generally below M Q but above Mf s for m ^ lGeV. The picture 
outlined above can be generalized to a continuous mass spectrum, and may also apply if fragmentation is 
replaced by formation of low mass 'satellites' in the course of a merger (i.e. if mergers produce, in addition 
to a large clump comprising most of the colliding mass, smaller ejected sub- lumps), although the value of 77 
could differ in that case. Needless to say, other pictures for the development of fluctuations might be more 
compelling than the one presented here, but may or may not lead to stationary statistics of the fluctuations. 

Although direct detection of elementary particle dark matter would be frustrated by dumpiness in 
the halo, indirect detection would be facilitated, perhaps to an embarrassing extent. Let us suppose that 
annihilation of dark matter particles produces some observable product - e.g. high energy gamma rays - with 
an efficiency Q(E)dE/E for energies near E. Then the (number) flux of observable annihilation products 
from a single clump of mass M a distance D from Earth is (using (av) ~ C/Mpim w 10 _27 (O /ig) _1 cm 3 s -1 ) 

dF cl (av)p cl (M)MC(E) 9 _ 2 _! S cl C(E)M/M Q 

E lLE 4^^> 3x10 cm & n fc§MGeV)m(pc)]' (3) 

where S c i — p c i /p, and the integrated flux from the halo of our Galaxy is 

dQ.dE Airm 2 O hl[m(GeV)] 2 v ' 

where (S c i) is the mean clump density contrast weighted by ip(M). (For predominantly p-wave annihilation, 
(av) oc (v 2 ) & 10~ 6 and Eqs. (3) and (4) would be correspondingly lower; the admixture of s- and p-wave 
in (av) is model-dependent.) For comparison, 5 the excess gamma ray flux above 100 McV measured by 
EGRET that is not correlated with Galactic HI column density is I ex = 1.47 ± 0.03 x 10~ 5 cm -2 s _1 sr _1 , 
and the isotropic component of this excess is estimated to be Ii SO = 1.5 ± 0.3 x 10~ 5 cm~ 2 s _1 sr _1 . Since the 
spectrum of the diffuse emission falls with increasing energy, 6 a halo flux as large as indicated by Eq. (4) 
could be excessive at photon energies ^ lGeV unless (5 c i)((E) is small. Also, EGRET has apparently found 
a small number of sources unidentified with any previously known astronomical objects. 7 These sources have 
total fluxes ~ 10 _7 cm _2 s _1 above 100 MeV, and no larger than ~ 10 _8 cm _2 s _1 above 1 GeV. Although 
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such fluxes could be consistent with Eq. (3) (depending on the value of 5 c i£(E) / D 2 ) it is probably unlikely 
that annihilation of multi-GeV particles would produce sub-GeV photons preferentially 8 and the mysterious 
sources could have a more mundane origin than annihilation in neutralino clumps. To the extent that 
small-scale concentrations of dark matter may be shown to be required physically, the absence of detectable 
annihilation products could be used to argue against the existence of certain neutralino candidates. 

This research was supported in part by NSF grants AST 91-19475 and AST 93-15375, and NASA grants 
NAG 5-3097 and NAG 5-2762. 
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